数値解析：理論的基盤と適切な定義

4次ルンゲ＝クッタ法やアダムス＝ムールトン公式といった数値ソルバーの力を発揮する前に、根本的な問いかけが必要です： 解は実際に存在するのか？また、安定しているのか？ 初期値問題（IVP）の理論的基盤は、数学的な「緑色信号」を提供し、離散化が意味のある物理的現実に収束することを保証します。これは、単なる数値ノイズに収束することを防ぎます。

基礎：リプシッツ連続性

誤差の伝播を制御するために、関数 $f(t, y)$ が「急激に変化しない」ことが必要です。この条件は リプシッツ条件によって形式的に定義されます。

定義 5.1：リプシッツ条件

関数 $f(t, y)$ が、集合 $D \subset \mathbb{R}^2$ 上で変数 $y$ に関してリプシッツ条件を満たすとは、ある定数 $L > 0$ が存在して、

$$|f(t, y_1) - f(t, y_2)| \le L|y_1 - y_2|$$

すべての $(t, y_1), (t, y_2) \in D$ に対して成り立つことを意味します。この定数 $L$ は、関数の縦方向の変化に対する「速度制限」となります。

例 1：リプシッツ定数の分析

集合 $D = \{(t, y) \mid 1 \le t \le 2, -3 \le y \le 4\}$ 上で関数 $f(t, y) = t|y|$ を考えます。平均値の定理（または絶対値の性質）により：

$|f(t, y_1) - f(t, y_2)| = |t|y_1| - t|y_2|| = |t| \cdot ||y_1| - |y_2|| \le |t| \cdot |y_1 - y_2|$.

領域内の $t$ の最大値が 2 であるため、リプシッツ定数は $L=2$ となります。

定義域の幾何学的整合性

穴だらけの定義域では初期値問題（IVP）を解くことはできません。我々は凸性によって形式的に定義されます。

定義 5.2：凸集合

集合 $D$ が凸であるとは、任意の二点 $(t_1, y_1)$ および $(t_2, y_2)$ に対して、以下の線分が $D$ 内に含まれることを意味します：

$$((1 - \lambda)t_1 + \lambda t_2, (1 - \lambda)y_1 + \lambda y_2)$$

ここで $\lambda \in [0, 1]$ です。これにより、解の経路の一部が有効な計算領域から「外れてしまう」ことがありません。

解の存在と一意性の定理

これらの条件が一致する場合、我々は 定理 5.4を適用します。関数 $f$ が凸集合 $D$ 上で連続であり、かつリプシッツ条件を満たすならば、初期値問題 $y' = f(t, y), y(a) = \alpha$ は 一意な 解 $y(t)$ を持つことになります。これは、オイラー法（$w_{i+1} = w_i + h f(t_i, w_i)$）のようなシンプルな手法から、予測補正ロジックまで正当化します：

$WC = w_{i-1} + \frac{h}{24}[9f(t_i, WP) + 19f(t_{i-1}, w_{i-1}) - 5f(t_{i-2}, w_{i-2}) + f(t_{i-3}, w_{i-3})]$.

🎯 核心原則：適切な定義

問題が 適切に定義された 一意な解が存在し、初期データに連続的に依存する場合を意味します。リプシッツ定数 $L$ が非常に大きい場合、問題は「剛性を持つ」となります。剛性方程式では、過渡成分は急速に減衰しますが、その導関数（大きさ $c^n e^{-ct}$）は減衰せず、そのため アルゴリズム 5.8：ニュートン反復付き台形則 安定性を維持するために必要です。

質問 1

定理 5.4 が初期値問題（IVP）の解の一意性を保証するために厳密に必要な条件はどれですか？

関数 $f$ は $D$ 全体で微分可能である必要がある

関数 $f$ は凸集合 $D$ 上で連続であり、かつリプシッツ条件を満たす必要がある

リプシッツ定数 $L$ は 1 より小さくなければならない

定義域 $D$ は非凸でなければならない

質問 2

初期値問題 $y' = 1 + t \sin(ty)$ （$0 \leq t \leq 2$）について、リプシッツ定数 $L$ の最良の推定値はどれですか？

$L = 1$

$L = 2$

$L = 4$

$L = 8$

質問 3

剛性方程式の文脈において、なぜ単純な陽的解法がしばしば失敗するのでしょうか？

過渡成分の減衰が遅すぎる

項自体よりも $n$ 階導関数の大きさ $c^n e^{-ct}$ が速く減衰せず、不安定性を引き起こす

凸性が破られている

剛性問題ではリプシッツ定数は常にゼロである

質問 4

ステップサイズ $h$ が $h = \sqrt{2\delta/M}$ よりも大幅に減少した場合、全体の誤差はどうなるか？

誤差は線形に減少し続ける

丸め誤差の支配が原因で、全体の誤差が増加する

解はより安定になる

リプシッツ定数は減少する

質問 5

『中点法』とはどのような方法ですか？

$w_{i+1} = w_i + h f(t_{i+1}, w_{i+1})$

テイラー法における $T^{(2)}$ を $f(t + h/2, y + h/2 f(t,y))$ で置き換える

アルゴリズム 5.8：ニュートン反復付き台形則

ニュートン後退差分公式の導関数

事例研究：数値的な適切な定義とダイナミクス

MATH007 ハイレベルチャレンジ

エンジニアは、質量 $m = 0.11$ kg の物体を鉛直上方に発射するモデルを作成しています。運動方程式は $mv' = -mg - kv|v|$ に従い、重力加速度 $g = 9.8$ m/s²、抵抗係数 $k = 0.002$ kg/m です。解く前に、数学的枠組みの妥当性を検証する必要があります。

課題 1

キルヒホッフの法則を安定性のアナロジーとして使用：$t=1.00$ から $1.04$ までの電流データ [3.10, 3.12, 3.15, 3.18, 3.24] が与えられたとき、中心差分法（$h=0.01$）を用いて、$L=2$, $R=10$ の条件下で $t=1.02$ における電圧 $\mathcal{E}(t) = L \frac{di}{dt} + Ri$ を近似します。

解答：
1. $di/dt \approx (i(1.03) - i(1.01)) / (2h) = (3.18 - 3.12) / 0.02 = 3.0$ を計算します。
2. $\mathcal{E}(1.02) = 2 \times 3.0 + 10 \times 3.15 = 6.0 + 31.5 = 37.5$ ボルト。

課題 2

初期値問題 $y' = y/t - (y/t)^2, y(1)=1$ （$[1, 4]$ 上）が適切に定義されていることを証明し、$t=2$ における実際の値を求めます。

解答：
1. $f(t,y) = y/t - (y/t)^2$. $\partial f/\partial y = 1/t - 2y/t^2$. $t \ge 1$ の有界領域では、この導関数は有界であり、定理 5.3 を満たします。
2. 実際の解は $y(t) = t/(1 + \ln t)$ です。
3. $y(2) = 2/(1 + \ln 2) \approx 2/1.6931 \approx 1.1812$。

課題 3

なぜ $u_2' = -24u_1 - 51u_2 - 9 \cos t + \frac{1}{3} \sin t$ といった剛性システムが、オイラー法（$w_{i+1} = w_i + h f_i$）のような陽的解法を実用不可能にするのか説明してください。

解答：
このシステムでは、ヤコビ行列の固有値の実部が非常に負です。過渡解は $e^{-ct}$ で $c$ が大きいほど急速に減衰します。陽的オイラー法が安定するためには、$h < 2/|\lambda_{\text{max}}|$ が必要です。ここでは $\lambda \approx -51$ なので、$h$ は極めて小さく（約 0.03）ならなければならず、$t$ が大きい場合の計算は非効率になります。代わりに アルゴリズム 5.8：ニュートン反復付き台形則 を使用すべきです。